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Abstract. We investigate anisotropic XXZ Heisenberg spin-1/2 chains with control fields acting on one 
of the end spins, with the aim of exploring local quantum control in arrays of interacting qubits. In this 
work, which uses a recent Lie-algebraic result on the local controllability of spin chains with "always- 
on" interactions, we determine piecewise-constant control pulses corresponding to optimal fidelities for 
quantum gates such as spin-flip (NOT), controlled-NOT (CNOT), and square-root-of-SWAP (VSWAP). 
We find the minimal times for realizing different gates depending on the anisotropy parameter A of the 
model, showing that the shortest among these gate times are achieved for particular values of A larger 
than unity. To study the influence of possible imperfections in anticipated experimental realizations of 
qubit arrays, we analyze the robustness of the obtained results for the gate fidelities to random variations 
in the control-field amplitudes and finite rise time of the pulses. Finally, we discuss the implications of our 
study for superconducting charge-qubit arrays. 

PACS. 03. 67. Ac Quantum algorithms and protocols - 03.67.Lx Quantum computation 
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1 Introduction 



Coherent control of quantum systems is one of the prereq- 
uisites for quantum information processing. While already 
simple arguments lead to the conclusion that almost any 
coupled quantum system can be controlled in principle [T] , 
the mathematical foundations of the subject are based on 
the notion of controllability and formulated using the lan- 
guage of Lie algebras [2J. In particular, a system is com- 
pletely controllable if its internal dynamics governed by 
external fields can give rise to an arbitrary unitary trans- 
formation in the Hilbert space of the system [3] . Both state 
control and the more general operator control have been 
implemented in a variety of systems [4] . 

Recent quantum control studies have focused their at- 
tention on interacting systems. A familiar example is fur- 
nished by spin chains, systems that can be used as data 
buses [5] for state- [51I7JI5JH] and entanglement transfer [TU] . 
In such systems, "always-on" interactions between the con- 
stituents (typically nearest neighbors) allow for a global 
control of the system dynamics by manipulating only a 
small subsystem, in the extreme case a single spin. The 
main question is then what is the smallest possible sub- 
system of a given system that one needs to act upon to en- 
sure the complete controllability, or, at least, the ability to 
perform certain pre-determined unitary transformations. 
This is the central idea behind the local-control approach. 



The fact that the local-control approach can be ad- 
vantageous in interacting systems provides an incentive 
for identifying minimal controlling resources that guaran- 
tee controllability in particular classes of systems. Quite 
recently, several Lie-algebraic results pertaining to local 
control of spin chains have been obtained [TTl[T2l fT3 , 14 , 
115] . For example, it was demonstrated that acting only on 
one of the end spins of an XXZ-Heisenberg spin chain 
ensures complete controllability of the chain [12] . Adopt- 
ing the last result as our point of departure, in this paper 
we investigate the feasibility of local operator control in 
qubit arrays modeled as spin-1/2 chains with Heisenberg 
interaction. In contrast to our recent proof-of-principle 
study |16j . where only the isotropic Heisenberg-coupling 
case was addressed, in the present work we discuss the 
case of (anisotropic) XXZ coupling. The main motiva- 
tion stems from the relevance of the XXZ-ca.se for imple- 
mentations of Josephson-j unction based superconducting 
qubit arrays [1711151 115] . 

We determine piecewise-constant control fields, acting 
only on the first spin in the chain, which lead to the highest 
possible fidelities for a selected set of quantum logic oper- 
ations: the spin-flip (NOT) of the last spin in the chain, as 
well as the controlled-NOT (CNOT) and the square-root- 
of-SWAP (VSWAP) gates applied to the last two spins. 
We optimize the gate fidelities with respect to the control- 
field amplitudes for three-spin chains. We then carry out a 
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sensitivity analysis, i.e., discuss the robustness of the ob- 
tained results with respect to random errors in the control 
fields, as well as finite rise/decay-times for control-field 
amplitudes. The present work is concerned with XXZ- 
Heisenberg spin chains and our conclusions apply to any 
physical realization of qubit arrays with this type of cou- 

pimg naag. 

2 System and method 

The total Hamiltonian of a Heisenberg spin- 1/2 chain of 
length N s reads 

H(t) =H + H c (t) , (1) 

where 

N s -1 

Hq = J {Si >X Si+l t x + Si t ySi+l t y + ASi :Z Si + l tZ ) , 

i=l 

(2) 

is a XXZ Heisenberg part with anisotropy A, and 

H c (t) = h x (t)S lx + h v (t)S ly (3) 

a Zeeman-like control part, with control fields h x (t), h y (t) 
acting only on the first spin. In what follows, we will also 
employ control Hamiltonians with fields in the x and z- 
directions. Whether the XXZ spin chain under consider- 
ation is ferromagnetic or antiferromagnetic is not crucial 
here, as we are concerned with operator control; aspects 
such as, for example, the different nature of the ground 
states in the two cases (separable vs. entangled) would 
only be consequential for issues related to, e.g., state con- 
trol or entanglement transfer. For definitcness, we will as- 
sume that J > and A > 0. It is useful to recall that 
the one-dimensional XXZ model has an antiferromag- 
netic ground state for A > 1, a ferromagnetic one for 
A < — 1, while for the intermediate values of A it is char- 
acterized by a critical gapless (quasi-long-range ordered) 
phase [2D]. 

For convenience, we hereafter set h = 1 and, in ad- 
dition, express all frequencies and control fields in units 
of the coupling strength J. Consequently, all times in the 
problem are expressed in units of 1/J. 

Since implementing control fields with a complex time 
dependence is difficult, we resort to piecewise-constant 
ones according to the following scheme. At t = we 
start acting on the first spin of the chain with an x con- 
trol pulse of amplitude h Xj i, which is kept constant until 
t = T. Thus the system is governed by the Hamiltonian 
= H + h Xt iS\ x . We then apply a y pulse with the 
amplitude h Vl i (Hamiltonian H y> i = Hq + h y< iSi y ) over 
the next interval of length T, etc. This sequence repeats 
until Nt pulses are carried out at t = tf = N t T. The full 
time evolution is described by 

U(tf) = Uy :Nt/2 U XtNt /1 . . . Uy A U X ,l , (4) 

where U x .i = e~' lH *- iT and U y .i = e~ lHy - iT are the respec- 
tive time-evolution operators corresponding to iT x ,i and 
H y _i, which can be evaluated using their spectral form. 



Our control objectives (target unitary operations) are 
both one-qubit gates, such as the spin-flip (NOT) opera- 
tion on the last spin of the chain Xn s := l<g>l<g>...(g>l(g)-X" 
(X being the Pauli matrix), and some entangling two- 
qubit gates. For instance, CN0Tjv 3 := 1 ® 1 ® . . . ® 
1 (£> CNOT performs the controlled-NOT operatio n on 
the last two qubits in the chain. Similarly, \/SWAPjv s := 
1 ® 1 ® . . . ® 1 ® VSWAP performs the VSWAP operation 
on the same pair of qubits. 

Unlike in many other control studies |11) , which make 
use of single-excitation subspaccs, we retain the full Hilbcrt 
space of the system. This puts constraints on the system 
size that can be treated within our framework. In what 
follows, we discuss three-spin chains. 

3 Controllability and reachability 

In Ref. |12) a very general graph-infection criterion was 
proven, which - as a special case - guarantees the complete 
controllability of XXZ Heisenberg spin chains through 
acting on one end spin. The more conventional approaches 
for proving complete controllability entail finding the di- 
mension of the relevant dynamical Lie algebra, a task for 
which special algorithms have been developed [21]. In the 
present problem, such an algebra is generated by the skew- 
Hermitian traceless operators {— iHq, — iS± x , — iS\ y } and 
has dimension d 2 — 1, where d = 2 Ns is the dimension of 
the Hilbert space of the system. Being generated by trace- 
less operators, this algebra is then isomorphic to su(d), 
the Lie algebra associated with the special unitary group 
SU(d) [2J]. 

Setting aside the issue of complete controllability, one 
might be interested to know if some particular unitary 
operations - on an otherwise not completely controllable 
system [23 , 2|j - are possible with an even smaller degree 
of manipulation, e.g.., a control field only in one direction. 
For such operations, equation (Q| goes over into U(tf) = 

U X} n ± ■ ■ ■ U Xt i. For example, the X^ s and v SWAP n b gates 
require only a control field in the x direction. To demon- 
strate this for Xn s , let C x be the dynamical Lie algebra 
generated by —iHq and —iS\ x , a subalgebra of su(d) with 
dimension 30 in a three-spin XXZ chain (note that the 
counterpart of this algebra in the isotropic-coupling case 
has smaller dimension, namely 18). For showing that Xm s 
belongs to the connected Lie subgroup e £x of SU(d) it 
suffices to find an element A G C x such that Xjy B = e A . 
Using the repeated commutators of the generators of C x , 
it can be demonstrated that X^ s is an element of this al- 
gebra. Xn 3 is both unitary and Hcrmitian, implying that 
Xjj = 1. It is then easy to show that A — —i^X^^ 
an element of C x , fulfills e A = —iXm s - Therefore, Xn, is 
r eachable using only an x control field. Recalling that the 
VSWAP gate on two qubits is given by [2"5] 

TSWAP = e*t e -*f (x®x+y®y+z®z) ^ (5) 

the reachability of a/SWAPjv, using an x control only 
readily follows from the fact that 1 ® 1 ® • • • <S> 1 S3 (X €5 
X + Y ® Y + Z ® Z) is an element of C x . 
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Fig. 1. Optimal control sequences for A = 5 realizing the (a) 
X z and (b) CNOT 3 gates with fidelity higher than 0.999. 



Table 1. Minimal times (in units of J ) needed to reach 
fidelities higher than 0.999 for the relevant gates in the x-y 
control case. The corresponding values in the x-z control case 
are given in the brackets. 



X, 



CNOT3 



1/SWAP3 



0.1 


23.3 


(19.0) 


73.2 


(60.2) 


14.2 


(14.2) 


0.2 


23.9 


(23.8) 


33.5 


(28.7) 


20.7 


(21.2) 


0.7 


18.2 


(18.0) 


22.2 


(22.6) 


18.9 


(16.8) 


0.8 


16.5 


(16.9) 


22.8 


(21.6) 


16.6 


(16.9) 


0.9 


15.4 


(13.9) 


18.4 


(18.2) 


1.6 


(1.6) 




1 K A 

lo.O 


(14. 9J 


1 T Q 

1 1 .3 


(16. ( ) 


1 r 
1 . 


(1.5) 


1.1 


15.4 


(15.8) 


17.3 


(16.4) 


1.5 


(1.5) 


1.2 


14.8 


(14.8) 


20.6 


(16.7) 


14.9 


(12.2) 


1.3 


16.2 


(16.0) 


20.0 


(19.0) 


15.2 


(14.8) 


2.0 


15.8 


(15.1) 


12.3 


(12.4) 


12.8 


(9.9) 


3.0 


12.0 


(13.4) 


12.0 


(12.2) 


10.2 


(8.8) 


4.0 


12.2 


(13.0) 


12.2 


(12.2) 


8.0 


(6.6) 


5.0 


12.2 


(13.2) 


11.2 


(11.1) 


5.4 


(5.3) 


6.0 


12.0 


(13.5) 


11.5 


(11.3) 


4.4 


(4.4) 


7.0 


12.9 


(13.9) 


11.8 


(9.8) 


3.8 


(2.0) 


8.0 


12.9 


(14.5) 


13.3 


(11.7) 


3.3 


(3.3) 


9.0 


13.6 


(15.0) 


13.1 


(10.4) 


4.3 


(3.0) 


10.0 


14.8 


(16.6) 


14.4 


(11.9) 


3.9 


(3.9) 


11.0 


21.7 


(17.9) 


13.1 


(10.8) 


2.4 


(2.4) 


12.0 


22.1 


(20.4) 


19.3 


(12.0) 


2.2 


(3.3) 



4 Target gates and minimal gate times 

In this section our goal is to find control fields leading to 
optimal fidelities for a chosen set of quantum gates, with a 
particular emphasis on minimal times needed for realizing 
different gates depending on the anisotropy A. 

In quantum operator control, the figure of merit is the 
gate fidelity 

F(t f ) = ±\tr[u\t f )U talgct ]\, (6) 

where U(tf) is the time-evolution operator of the sys- 
tem at time t — tf (Eq. Qj) and /7t arg ct stands for the 
quantum gate that we want to realize. We perform opti- 
mization, i.e., maximize the gate fidelity with respect to 
the N t control-field amplitudes, for varying number (N t ) 
and durations (T) of pulses (hence different total evolu- 
tion times tf). We make use of a quasi- Newton method 
due to Broyden, Fletcher, Goldfarb, and Shanno (BFGS- 
algorithm) [2B]. It should be stressed that, much like other 
optimization approaches, this algorithm ensures only con- 
vergence to a local maximum. Therefore, to determine a 
globally-optimal sequence of control-field amplitudes (for 
a given target gate and given value of A) we ought to re- 
peat the optimization process for a number of different ini- 
tial guesses for these amplitudes. We generate these initial 
guesses using a uniform random number generator [26 . 

An alternative to fixing the pulse durations and maxi- 
mizing over the control-field amplitudes would be to keep 
the amplitudes constant and treat the pulse durations as 
variable control parameters. However, we choose optimiza- 
tion over the control-field amplitudes since this approach 



allows us to easily fix tf and determine its minimal value 
for implementing the desired gate for any fixed value of 
the parameter A. 

The obtained results for the gate fidelities have the fol- 
lowing two salient features. Firstly, for fixed parameters of 
the model and fixed total evolution time t / , the fidelity for 
any given gate can increase significantly with increasing 
N t (or, equivalently, decreasing T). In other words, more 
rapid switching leads to higher fidelities. For instance, in 
the case of the CNOT 3 gate with A = 1.3 and tf = 30, 
for N t = 10, 20, 30, 40, 50, 60, 70 we obtain the respective 
fidelities F = 0.455, 0.697, 0.837, 0.953, 0.995, 1 - 10~ 4 , 1 - 
10~ 8 . Secondly, for each gate there exists a minimal value 
of i/ (i.e., minimal gate time), below which fidelities close 
to unity cannot be reached regardless of the value of Nt ■ 
The obtained minimal gate times for different values of A 
in the x-y (x-z) control cases are given in Table [1] Ap- 
parently, there exists an optimal value of A which cor- 
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responds to the shortest among these times. For the A3 
and CNOT3 gates, for example, these values are around 
A = 5. The corresponding optimal sequences of x and y 
control pulses for the X 3 and CNOT3 gates are shown 
in Figures [ITa) andQJb), respectively. Since ideal steplike 
pulses cannot be realized in practice, in Ref. [TB] we also 
studied frequency-filtered control fields and showed that 
sufficiently high fidelities can still be retained. 



In Ta ble [2] the minimal times are given for the A3 and 
VSWAP 3 gates realized using only an x control field. It is 
interesting to compare these minimal times to the above 
case with both x and y (or x and z) controls. For small 
values of A (with the exception of A = 1) the minimal 
times for realizing the X3 gate in the x-only control case 
are significantly longer than their counterparts in the x-y 
(x-z) case. In contrast, for larger A these times become 
more and more similar. Finally, for A > 11 the minimal 
times in the x-only control case are even shorter than in 
the x-y and x-z cases. Since £-only control is easier to 
implement, this surprising observation provides an addi- 
tional argument for using x-only control in the regime of 
interest for superconducting charge qubits. 

As is well known [35] , the vSWAP gate on two qubits 
is naturally implemented by the isotropic Heisenberg Hamil- 
tonian after a time r = 7r/2 w 1.57. As can be seen in 
Tabic [TJ the minimal voWAP 3 -gate times indeed seem to 
correspond to A ~ 1 and are only slightly shorter than 
in the control-free case. This is despite the fact that our 



Table 2. Minimal times (in units of J -1 ) needed to reach 
fidelities higher than 0.999 for the relevant gates using only an 
x control field. 
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VSWAP3 


0.1 
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80.6 


3.0 


19.2 
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0.2 
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45.0 
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18.4 


9.5 
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25.0 
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37.1 
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Fig. 2. (Color online) Average fidelity versus half-width (5) for 
the v / SWAP 3 gate with A = 1.2 and fixed total time t f = 60. 



%/SWAP 3 gate performs the VSWAP operation on the last 
two qubits (leaving the state of the first qubit unchanged) 
while the Heisenberg Hamiltonian of equation @ also con- 
tains the interaction between the first two qubits. Thus we 
can conclude that the role of control fields in this case is 
to counteract the effect of the free evolution of the first 
qubit governed by Hq . 

Generally speaking, the minimal gate times can in prin- 
ciple be found based on the time-optimal unitary opera- 
tion formalism put forward by Carlini and co-workers [271 ■ 
This method requires solving a system of coupled nonlin- 
ear equations for Lagrange multipliers resulting from the 
quantum brachistochrone equation. In practice, extracting 
minimal times for different quantum gates in this way is 
feasible only when the time evolution of the total Hamilto- 
nian of the system is as simple as to allow for an analytical 
solution of these equations. This is possible, for instance, 
when this Hamiltonian has a block-diagonal form in the 
computational basis, where each block commutes with it- 
self at different times. In the problem at hand this is not 
the case, therefore an alternative strategy for finding min- 
imal times is required. 



5 Robustness to random errors and finite 
pulse rise times 



In the following, we analyze the sensitivity of the fidelity 
to random errors in the control-field amplitudes, as well 
as to a finite rise time. 

The random errors in control-field amplitudes are as- 
sumed to follow a uniform distribution of half- width 8. For 
given 8, we generate a large sample of N ~ 1000 control 
fields affected by random noise, for which we recalculate 
the fidelity. We are interested in the behavior of the av- 
erage fidelity F = Yli=i Fi/N, where the Fi are fidelities 
for specific realizations of the random field, versus 6 for 
the gates of interest and varying values of A. 
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In our previous work |16j . using the isotropic Heisen- 
berg model (A — 1) as an example, it was demonstrated 
that the shape of the fidelity decay curves (F vs. S) de- 
pends on the number of control pulses N t and their length 
T. Provided that the system satisfies the conditions for 
complete controllability, the saturation regime of the av- 
erage fidelity sets in for S > 2J. The universal saturation 
value is l/<i, where d is the dimension of the Hilbert space 
of the system. Importantly, for fixed tf = N t T, the av- 
erage fidelity is closer to the intrinsic (in the absence of 
random errors) optimal values for larger N t (faster switch- 
ing), this being a consequence of general properties of 
systems that exhibit competition between the resonance- 
and relaxation- type behavior [IB]. Therefore, more rapid 
switching leads not only to higher intrinsic fidelities in 
the absence of randomness (recall section 2]), but also ren- 
ders these fidelities less sensitive to random errors. This is 
a manifestation of an intrinsic robustness of the system. 
Figure [3] illustrates that these features are also present in 
the anisotropic XXZ case. 

The sensitivity to random errors in the control-field 
amplitudes depending on the anisotropy A is illustrated 
in Figure [3] As can be inferred from this figure, for larger 
A the system is less sensitive to random errors. 

Another unavoidable source of imperfections in qubit- 
array realizations is the finite rise time of the control fields. 
Instead of a stepwise behavior, experimental control fields 
hj, n (J — x,y; n = 1, . . . , N t /2) are expected to have a 
finite rise/decay time t. Figure [4] shows the dependence of 
the fidelity on the finite rise time. For larger values of A, 
the fidelities of optimal control sequences seem to be less 
affected by the finite rise time. 

The central result of this section is that values of A > 1 
lead to both shorter gate times and a reduced sensitivity 
of the fidelity to random errors in control field and finite 
rise time. This provides a guiding principle for future im- 
plementations of qubit arrays. 



6 Discussion and Conclusions 

Our results are of direct relevance to superconducting 
qubit arrays [BJ. One-dimensional Josephson arrays of ca- 
pacitively coupled superconducting islands can be described 
as XXZ Heisenberg spin-1/2 chains [28,29]. In general, 
the A Y"-part of Hamiltonian is characterized by a nearest- 
neighbor interaction, whereas the Z-part will also have 
coupling contributions beyond nearest neighbors. How- 
ever, by properly choosing the junction capacitances and 
the capacitance of each island to the back gate of the struc- 
ture, the Z-part will also be approximately of nearest- 
neighbor type. The correspondence between the parame- 
ters of the Josephson array and the spin chain is as fol- 
lows: the Josephson energy Ej of the junctions coupling 
the islands corresponds to the exchange coupling constant 
J of the spin system and can be controlled by a magnetic 
field if we assume that the coupling junctions are realized 
as SQUIDs. The parameter AJ of the spin system corre- 
sponds to the charging energy Eq, i.e., the anisotropy pa- 
rameter A corresponds to Eq/Ej. Values of A like those 
studied in Tables [1] and [2] can be experimentally realized. 
Finally, the first island should form a charge qubit, and 
the control field h z corresponds to the gate voltage, while 
h x and h y play the role of the Josephson energy of this 
charge qubit. Our study shows that, in principle, arbitrary 
quantum algorithms can be realized on a one-dimensional 
Josephson array by controlling only the first island in the 
array. 

In summary, we have shown that local control of the 
first spin of an anisotropic XXZ Heisenberg spin-1/2 chain 
enables universal quantum computation. Using a recent 
Lie-algebraic result on the local controllability of spin chains 
with "always-on" interactions, we have determined con- 
trol pulses leading to optimal fidelities for quantum gates 
such as spin-flip (NOT), controlled-NOT (CNOT), and 
square-root-of-SWAP (VSWAP). We have found the min- 
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N t = 70 and T = 1 corresponding to the CNOT3 gate. respond to Nt = 70 and T = 1. 
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imal times for realizing different gates depending on the 
anisotropy parameter A of the model, showing that the 
shortest among these gate times are achieved for particu- 
lar values of A larger than unity. Another surprising re- 
sult was that in the regime of interest for superconduct- 
ing charge qubits, the minimal times in the simpler x-only 
control case can be even shorter than in the x-y and x- 
z control cases. We have also analyzed the sensitivity of 
the obtained results for the gate fidelities to random vari- 
ations in the control-field amplitudes and finite rise time 
of the pulses. Our results are independent of a partic- 
ular experimental realization of the XXZ chain, yet, a 
superconducting Josephson array would be a particularly 
appealing candidate. Our investigation paves the way for 
future studies, involving more sophisticated control strate- 

gies pnusu- 
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